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The dynamics of instantons and anti-instantons in lattice QCD can be studied by analysing the action and 
topological charge of configurations as they approach a self-dnal or anti-self-dual state, i.e. a state in which 
S/So = IQ I- We use cooling to reveal the semi-classical structure of the configurations we study. Improved 
actions which eliminate discretization errors up to and including 0{a^) are used to stabilise instantons as we cool 
for several thousand sweeps. An analogously improved lattice version of the continuum field-strength tensor is 
used to construct a topological charge free from 0{a^) discretization errors. Values of the action and topological 
charge obtained with these improved operators approach mutually-consistent integer values to within a few parts 
in 10"^ after several hundred cooling sweeps. Analysis of configurations with |Q| « 1 and |Q| « 2 supports the 
hypothesis that a self-dual |Q| = 1 configuration cannot exist on the 4-torus. 


1. INSTANTONS 

In 1975 Belavin et al. 0 determined that the 
vacuum of QCD contained non-trivial local min¬ 
ima of the Euclidean action, 

Se = ^ j ( 1 ) 

Since this action is required to be finite in the case 
of semi-classical solutions (i.e. solutions obtained 
by minimizing the action) the definition of the 
field-strength tensor 

Ffiv = — dyAfj^ + [^AD (2) 

implies that 

^ U-^d^U (3) 

as [a;I goes to infinity. The gauge transformation 
U is interpreted as a mapping from —> SU (2), 
that is, a mapping from the sphere at spacial in¬ 
finity to the gauge group SU{2) of the transfor¬ 
mation U {SU (2) results can be readily extended 
to SU (3)). The mapping U can then be shown to 
have a winding number given by the equation 

n = j d/^xFvF^yF^y (4) 

and hence we can associate a winding number 
with each vacuum state. Instantons are inter¬ 
preted as tunnelling amplitudes between vacuum 


states with different winding numbers, and their 
topological charge is given by the difference be¬ 
tween the winding numbers of the states that 
they connect. Further important properties of 
(anti-)instantons are spherical symmetry in the 
four dimensions of Euclidean spacetime, topolog¬ 
ical charge of ±1, and an associated action of 
So = which is independent of the radius 

of the instanton. 

2. CONNECTING THE LATTICE TO 
THE CONTINUUM 

The lattice is a 4-dimensional approximation 
to Euclidean spacetime. This means that physics 
on the lattice should ideally approach physics in 
the continuum if we can take the limit of infinite 
lattice volume and infinitessimal lattice spacing 
{V —> oo and a —> 0). Unfortunately reducing 
the lattice spacing indefinitely is not practical, 
nor is increasing the total lattice volume, since 
these approaches to the continuum limit rapidly 
increase the time taken to perform computations. 
We have chosen to remove the non-physical ‘edge’ 
of the lattice (simulating an infinite volume) by 
defining periodic untwisted boundary conditions 
on our lattices. We are therefore performing 
numerical QCD on a 4-toroidal mesh of points. 
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The discretization errors introduced by the fi¬ 
nite lattice spacing are reduced algebraically by 
Symanzic improvement and non-classical er¬ 
rors arising from the self-couplings of the gluon 
fields are dealt with by the use of mean-field im¬ 
provement 1 ^. 

Since the lattice we use is topologically equivalent 
to a 4-torus, we may consider the consequences 
of the Nahm transform [Q for the stability of 
the configurations we study. The Nahm trans¬ 
form is a duality mapping which interchanges be¬ 
tween an SU{N) configuration with topological 
charge Q on the torus and an SU{Q) configura¬ 
tion with topological charge N on the dual torus. 
Since there are no instanton solutions in C/(l), the 
Nahm transform implies Q||] that there can be 
no self-dual IQI = 1 configurations on the torus. 


3. THE LATTICE ACTION AND TOPO¬ 
LOGICAL CHARGE 


3.1. The Wilson Action 

The Yang-Mills action on the lattice Q is 


Swil = /3 EE 1 

X ^<iy 


1 (ReTrW^ixi)(x)) 


( 5 ) 


where is the m x n Wilson loop in the 

pL — V plane. When we expand the plaquette 
around the point x and set j3 = &lg'^ it 
can be easily shown that (||) approaches the con¬ 
tinuum Yang-Mills action, plus 0{a?) corrections 

5wn + ©(a^) (6) 

in the continuum limit a —> 0 . 


• L( 2 . 2 ) = 

• L(b3) = 1 

• L(3.3) = 

A general improved action can then be written in 
the form 

C4SiL^^'^^)+C5S{L^^’^^). (7) 

where the ci,..., C 5 are improvement factors which 
take the values 

Cl = (19-55 c5)/9 

C2 = (1-64c5)/9 

C 3 = (640c5 - 64)/45 
C 4 = 1/5 — 2c5 

and C 5 is a free parameter which we can use to 
“tune” the action. By setting C 5 = 0 we cre¬ 
ate a 4-loop improved action, denoted 4LIS for 
short, and by setting C 5 = 1/10 we eliminate the 
contribution of and thereby 

creating a 3-loop improved action, denoted 3LIS 
for short. For the purposes of investigating a 5LIS 
we choose to follow DeForcrand et al. and set 
C 5 = 1/20, midway between the 3LIS and 4LIS 
values. It is important to note that the 3LIS and 
4LIS are just special cases of the general 5LIS, 
and all choices of C 5 are free from discretization 
errors up to and including 0{a'^), however their 
0{a^) errors will be different, leading to slightly 
different results. 


3.2. Improving the Action 

We can eliminate 0{a^) errors by constructing 
the Wilson action from a linear combination of 
the plaquette and the average of the 1 x 2 and 
2x1 rectangular Wilson loops. We may of course 
choose to employ other loops to eliminate higher- 
order error terms. DeForcrand et al. have 
previously constructed an action free from G(a^) 
and 0{a‘^) errors by using the five planar Wilson 
loops 

• l(i4) = 


3 . 3 . Topological Charge 

On the lattice the topological charge is calcu¬ 
lated as the sum of the local topological charge 
density over all lattice sites, 

Q — ^ ^ 9(2^) = 32 ^ ^ ^ C/ii/po-TrF/jyEpCT- ( 8 ) 

X X 

If we have a dilute instanton gas with n/ instan- 
tons and anti-instantons the total topological 
charge of the configuration will be 

Q = m - riA- 


( 9 ) 
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We hence desire the topological charge calculated 
on the lattice to approach integer values as the 
configuration being studied is cooled towards a 
dilute instanton gas condition. Unfortunately dis¬ 
cretization errors tend to lead to non-integer val¬ 
ues of Q. Fortunately, however, we can improve 
the topological charge. 


To construct an improved field strength tensor we 
utilise the same five planar Wilson loops that are 
used in the construction of the improved action, 
but in this case we calculate from the 

clover average of four m x n Wilson loops in the 
^ — V plane. The improved field strength tensor 
can then be written 


3.4. The Improved Field-Strength Tensor 

To improve the topological charge we have cho¬ 
sen to improve the field-strength tensor di¬ 
rectly, by analogy with the improvement of the 
action, and substitute this into Eq. (||), the defi¬ 
nition of the topological charge. 

Consider the full expansion of the plaquette 

= e^3§Adx 

= 1 + ig ^ Adx — ^ (j) AdxY 

+0{g^) 

= l+ig + 0{a'^)] - 

+0(a^g3). (10) 

We wish to extract the second term (within the 
square parentheses) on the last line. This may be 
achieved by making the following construction; 


= l + ^gIAdx-^{<j 

[ Adxf 


+0{g% 

= l-igfAdx-^{^ 

^ AdxY 


+0{g^). 


( 11 ) 

Writing 

A = IU*^^xi) _ 

fii/ fiy 


( 12 ) 


we may obtain a term proportional to the field 
strength 


A — -Tr^^ = g Adx + 0{g^) 

= ga^Ff,^ P 0{ga^). 

(13) 

Notice that in order to enforce the tracelessness of 
the Gell-Mann matrices we have subtracted one- 
third of the trace of A. 


+k5F^.{L^^’^^) (14) 

where the ki,...,k 5 are improvement factors 
which take the values 

fci = 19/9-55^5 
^2 = 1/36 — 16^5 

k3 = 64^5 - 32/45 
A:4 = 1/15 — 6^5 

and fca is a free parameter, which enables us to 
create 3-loop, 4-loop, and 5-loop improved defini¬ 
tions of the field-strength, all of which are 0{a‘^)- 
improved. 

In addition to constructing Q, the field strength 
may be used in (||) to create a “reconstructed ac¬ 
tion” . Since this action is improved in a different 
manner to the cooling action, comparison of both 
actions may serve as a mechanism with which to 
test the improvement procedures. 

In this brief report we will simply “test the wa¬ 
ter” with the values 05 = ^5 = 1/20 (5-loop), C 5 = 
^5 = 0 (4-loop), and C 5 = 1/10, k^ = 1/90 (3- 
loop), and show that we obtain excellent results 
even though many other 5-loop operators, corre¬ 
sponding to other choices of C 5 and fcs, are possi¬ 
ble. 

4. RESULTS AND COMPARISON 

4.1. Constructing Configurations 

The configurations used in our investiga¬ 
tions are generated using the Cabibbo-Marinari 
pseudo-heatbath algorithm |^. We generate 
SU (3) configurations by looping twice over three 
diagonal SU{2) subgroups. From a cold start (all 
link values set to zero) we thermalize for 5000 
sweeps using an (!l(a^)-improved action and a 
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Figure 1. Action obtained with 1-loop (a) and 
2-loop (o) cooling. 


fixed mean-link value, and then select configura¬ 
tions every 500 sweeps to create an ensemble of 
configurations with non-trivial topology, which 
are statistically distinct from each other. 

To reveal the semi-classical structure of the con¬ 
figurations we use a cooling algorithm which 
likewise uses diagonally embedded SU{2) sub¬ 
groups, with appropriate link partitioning |l^ . 
The results we present are produced from con¬ 
figurations on a 12^ X 24 lattice with f3 = 4.60, 
providing a lattice spacing of a = 0.125 fm. 

4.2. Comparison of Cooling Schemes 

In practice, because of the large discretization 
errors present, a cooling scheme based upon the 
plaquette action will destroy topological struc¬ 
ture over the whole lattice if it proceeds for long 
enough. However, a cooling scheme using an 
improved action will stabilise at a state with 
non-trivial structure. We would anticipate cool¬ 
ing schemes with small discretization errors will 
reproduce expected continuum behaviour, espe¬ 
cially integer values of Q and S/Sq, and remain 
stable for hundreds of cooling sweeps. When com¬ 
paring several cooling schemes, the one which 
fulfills these criteria best will be the preferred 
scheme for studying gauge field dynamics. 

In Figure we can see how the use of a 2- 
loop, O(a^)-improved action in the cooling algo¬ 
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Figure 2. Action obtained with 2-loop (a), 3LIS 
(a), 4LIS (marked), and 5LIS (•) cooling. 


rithm stabilises a configuration at a state with 
S/So ~ 6, indicating the presence of six in- 
stantons and/or anti-instantons, while plaquette- 
based cooling gradually destroys the structure of 
the configuration. 

In Figure ^ we compare the values of the action 
obtained with 2-loop, 3LIS, 4LIS, and 5LIS-based 
cooling on the same configuration used in Fig¬ 
ure 1^. We can see that the 4LIS curve drops well 
below the desired integer value of six, indicat¬ 
ing the presence of large negative discretization 
errors. This comparison and similar work with 
other configurations have lead us to conclude that 
3LIS-based cooling gives the most integer-like re¬ 
sults in general, and furthermore it is faster than 
4LIS or 5LIS-based cooling, since it requires fewer 
computational steps. 

4.3. Comparison of Topological Charge 
Improvement 

In Figure || we compare the various improved 
Q values obtained as we cool a configuration. 
The 3-loop improved operator gives the most 
continuum-like results, coming to within around 
1 part in 6000 of an integer, and has the small¬ 
est computational cost of all the C(a‘^)-improved 
operators. Similar comparisons with other config¬ 
urations have lead us to select the 3-loop operator 
as our preferred choice for studying the topologi¬ 
cal charge of cooled configurations. 
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Figure 3. \Q\ of 2L (x), 3L (a), 4L (□), and 5L 
(o), with 3LIS cooling 


4.4. Investigation of the Nahm Transform 

Having established our preferred operators for 
investigating instanton dynamics, we are now in 
a position to compare the stability of IQj = 1 and 
Ql ^ 1 configurations, to determine if we can find 
any evidence that self-dual \Q\ = 1 configurations 
are not permitted on the 4-torus, as implied by 
the Nahm transform. Figure ^ shows an over¬ 
lay of the reconstructed action and topological 
charge of CFG 1, a |(5| « 1 configuration, and the 
corresponding values for CFG 2, a |(5| ~ 2 con¬ 
figuration (shifted down by an increment of one 
unit so that they overlap with those of CFG 1 on 
the scale of this diagram). These configurations 
are separated by 8000 thermalization sweeps. We 
can clearly see that the |Q| ~ 2 configuration sta¬ 
bilises quickly at a self-dual near-integer value, 
while the other configuration does not achieve 
self-duality and eventually destabilises. Investi¬ 
gations of other configurations have shown equiv¬ 
alent behaviour. 

5. CONCLUSIONS 

Highly improved operators enable us to cool 
configurations to stable non-trivial self-duality. 
Counting the number of instantons and anti- 
instantons present by evaluating S/Sq and Q al¬ 
lows us to detect evidence of continuum-like be¬ 
haviour. The long-term stability of \Q\ « 2 con¬ 
figurations indicates that our cooling algorithms 
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Figure 4. 3-loop improved |Q| and reconstructed 
S/Sq of CFG 1 and GFG 2 


are highly reliable, and hence the destabilization 
of I <51 ~ 1 configurations is due to the conse¬ 
quences of the Nahm transform, indicating that 
on the lattice, just as in the continuum, a con¬ 
figuration may not simultaneously fulfill the two 
criteria S/Sq = 1 and \Q\ = 1 on the 4-torus. 
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